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Low temperature specific iieat has been measured for an overdoped La2-xSrxCuOi ( x = 0.22 ) 
single crystal. The quasiparticle density of states ( DOS ) in the mixed state is found to deviate from 
the predicted scaling law Cvoi ~ H f{T / H^^^). However, this scaling behavior is nicely reconciled 
if one considers the normal core region ( ^ ~ 2lJi ) which gives small contribution to the total 
DOS. The radius of the core and other parameters derived are consistent with reported values. Our 
results suggest that there is no zero-bias conductance peak ( ZBCP ), which is predicted by the 
simple Bogoliubov de-Gennes theory in the vortex core of a d-wave superconductor. 
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One of few points with consensus in the cuprate super- 
conductors is the dx2-y2 pairing symmetry in hole doped 
cuprates. This has been supported by tremendous exper- 
iments Q both from surface detectionjl H S 0, and 
bulk measurements[7ll5 liL llo| . In a d-wave superconduc- 
tor with line nodes in the gap function, the quasiparticle 
density of states ( DOS ) N{E) rises linearly with en- 
ergy at the Fermi level in zero field, N{E) cx — Ep\ , 
resulting in an electronic specific heat Cg = aT"^ |l2l |. 
where a = 2.527„/Tc and 7„ is the specific heat coef- 
ficient reflecting the DOS at the Fermi level of normal 
state. In the mixed state with the field higher than a 
certain value, the DOS near the Fermi surface becomes 
finite, therefore the quadratic term Ce — aT^ will be sur- 
passed and substituted by both the localized excitations 
inside the vortex core and the de-localized excitations 
outside the core. Volovik pointed out that for d- 
wave superconductors in the mixed state, supercurrents 
around a vortex core lead to a Doppler shift to the quasi- 
particle excitation spectrum, which affects strongly the 
low energy excitation around the nodes. It was shown 
that the contribution from the delocalized part will pre- 
vail over the core part and the specific heat is predicted 
to behave as C^oi — k^nT H / Hc2 with k a parameter 
in the order of unity. This prediction has been verified 
by many specific heat measurements which were taken as 
the evidence for d-wave pairing svmmetrvjsl 1^ HoL 
In the finite temperature and field region a scaling law is 
proposed^ 
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where / is a unknown scaling function. This scaling law 
has been proved in YBCO 9, 10] and in optimally doped 
LSC0[l2|. It remains however unclear whether this scal- 
ing law is still valid in the overdoped region. 



Another important but controversial issue is the vor- 
tex core state in the cuprate superconductors. By solv- 
ing the mean-field Bogoliubov-de Gennes ( BdG ) equa- 
tion, theoretically it is suggested that a zero-bias con- 
ductance peak ( ZBCP ) exists in the vortex core [l5lll6| . 
However this is in sharp contrast with the experimen- 
tal observations [H H E US IH i^. The absence 
of ZBCP was attributed to the presence of id^y or is 
componentsji^ . In this Letter we show that the DOS due 
to vortex quasi-particle excitations deviates from Simon- 
Lee scaling law, but it can be nicely reconciled if the 
vortex core size is taken into account. Furthermore a 
consequence from our analysis is the absence of ZBCP in 
the vortex core and the disappearance is intrinsic rather 
than due to the second component of order parameters. 

The single crystals measured in this work were pre- 
pared by the travelling solvent floating-zone technique. 
In this Letter we present data of an overdoped single 
crystal ( x — 0.22, Tc — 27 AK) as characterized by 
AC susceptibility and DC magnetization ( shown by the 
inset in Fig.l ). The quality of our sample has also 
been characterized by x-ray diffraction patterns show- 
ing only (001) peaks, and resistive measurement show- 
ing a rather narrow transition temperature ATc < 2 K. 
A piece about 28.55 mg in mass, 3.66 x 2.3 x O.Smm'^ 
in dimension, was chosen for the specific heat measure- 
ment. The heat capacity d ata p resented here were taken 
with the relaxation method [23 based on an Oxford cryo- 
genic system Maglab. The heat capacity is determined 
by a direct measurement of the thermal time constant, 
T — (C+Cadd)/ here C and Cadd are the heat capacity 
of the sample and addenda ( including a small sapphire 
substrate, small printed film heater, tiny Cernox temper- 
ature senser, 025 fim gold wire leads, Wakefield thermal 
conducting grease ( lOOfig ) ) respectively, where 
the thermal conductance between the chip and a thermal 
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FIG. 1: Specific heat coefficient C/T vs. at magnetic 
fields ranging from to 12 T for the overdoped single crystal. 
The inset shows the diamagnetic transition at around 27.4 K 
determined by the crossing point of the extrapolating line of 
the most steep part with the normal state background M = 
0. 

link. The value Cadd has been measured and subtracted 
from the total capacitance, thus C value reported here is 
only the capacitance of the sample. We have also checked 
the field dependence of Cadd and found that the change 
( if any ) of Cadd under 12 T is in the same order of the 
noise background here ( 20 nJ/K at 5 K and 40 nJ/K at 
20 K ). The influence of the magnetic field ( 12 T ) on 
the readout of the thermometer is below 0.02 K and can 
be neglected. During the measurement the field was ap- 
plied parallel to c-axis and the sample was cooled to the 
lowest temperature under a magnetic field (field-cooling) 
followed by data acquisition in the warming up process. 

Fig.l shows the specific heat coefficient C/T as a func- 
tion of in magnetic fields ranging from to 12 T. 
The separation between each field can be well deter- 
mined. In low temperature region the curves are rather 
linear showing that the major part is due to phonon con- 
tribution Cph = (3T^ ■ In addition the curve at zero 
field extrapolates to a finite value ( 70 ) at K in- 
stead of zero as observed in other cuprate supercon- 
ductors. This may be interpreted as potential scatter- 
ing to the vanishing gaps near the node of dx^-y'^ gap 
function due to small amount impurities As also ob- 
served by other groups for La — 214 system, the anoma- 
lous upturn of C/T due to the Schottky anomaly of free 
spins is very weakpr]. This avoids the complexity in the 
data analysis. It is known that the phonon part has a 
very weak field dependence, this allows to remove the 
phonon contribution by subtracting the C/T at a cer- 
tain field with that at zero field. The results after the 
substraction are shown in Fig. 2. The subtracted values 
A7 = 7H - 70 = C{T,H)/T - C{T,H = 0)/T show a 
rather linear T behavior in low temperature region ( be- 
low 8 K ) as indicated by the solid lines. Our analysis is 



FIG. 2: Temperature dependence of the subtracted specific 
coefficient 7^f - 70 = C/T{H) - C/T{H = 0). A linear be- 
havior is clearly seen in low temperature region which is not 
in accord with the proposed scaling law (eq.l) by Simon and 
Lee ( see text ). The linear lines in low temperature region 
are guides to the eyes. From these lines one can determine the 
zero temperature intercept A7 and the slope d'y/dT shown in 
Fig.3 

based on the data below 12 K. In the following we will 
show that the field dependent slope of the linear part in 
low temperature re gion shown in Fig. 2 directly deviates 
from the Simon-Lee [iJI scaling law ( eq.l ). 

According to the scaling law, the low temperature ex- 
pansion of right hand side reads 

Croi=baH + biTVH + b2T^ + o{T^) (2) 

where bo=0 because Cvoi/T should be finite when 
T — and H 0. Since o{T^) is very small in low tem- 
perature region, one has Cyoi/T = bi^/H + b2T. When 
Hci << H << Hc2, the total specific heat contains four 
parts: Doppler shift term from the region outside the 
core Cvoi, the inner vortex core term Ccore oc HT^ the 
small impurity scattering term 70T and the phonon term 
Cph- At zero field, the total specific heat contains three 
parts: the small impurity scattering term 70T and the 
phonon term Cph both depend on magnetic field weakly, 
and a quadratic term aT^ due to the thermal excitation 
near the nodal region. Thus A7 can be written as: 

A-f = -fH -10 = blVH +{b2-a)T + -faoreH (3) 

From eq.3 it is clear that A7 depends on T through the 
second term, however the slope 62 — a is independent 
according to the definition. This clearly indicates that 
the Simon-Lee scaling law cannot be directly applied to 
interpret the field dependent slope of A7 vs. T as shown 
in Fig.2. 

In order to understand the underlying physics, still 
based on the Simon-Lee scaling law, we propose that the 
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core size effect has a sizable influence on the total vortex 
quasi-particle excitations. By taking account the vortex 
core size (2^), one can rewrite A7 as: 

A7 - {hl^/H + h2T) X (1 - eVi?2) _ ^ ^^^^^jj 

where ^ is the radius of the normal core, Ra is the 
radius of a single vortex — (pQ/nH . Thus eq.4 can be 
written as: 

A7 = 6ix/i7x {l-^H) + {h2-a)T-h2^HT+-f,oreH 

(5) 

One immediately realizes that the third term in eq.5 is 
just what we need for interpreting the difficulty as men- 
tioned above. Next let us have a closer inspection at 
the data and derive some parameters. At zero tempera- 
ture, only the first term and the last term are left. The 
values of A7(r = 0) are determined from the extrapo- 
lation of the linear lines in Fig. 2 to K and presented 
in Fig. 3. The solid line is a fit to the data using the 
first term in eq.5 yielding hi = 1.9 ± 0.042m Jif-^T-^/^ 
and 7rCV0o = 0.0067 ± 0.002 and thus ^ = 2ll. The 
value ^ = 21A derived here is quite close to that found 
in Nernst jj^ and STM measurementsj^. We also tried 
to use the first term together with the last term to fit 
the data but find out that the contribution from the 
last term is extremely small. The first term here de- 
scribes the zero temperature data very well, indicating 
the absence of a bulk second order parameter such as 
idxy or is since otherwise the nodal point would be filled 
completely and the Doppler shift had very weak effect 
on the quasi-particle excitations. The inset of Fig. 3 
shows the field dependence of the slope of the linear 
part in Fig. 2. It is clear that the slope increases with 
H above 1 T. This can be exactly anticipated by the sec- 
ond and third terms in eq.5. From the inset of Fig. 3 
one obtains a — 62 = 0.21m Jmol^^ K^^ and b2Tr£,'^ / (j)o ~ 
0.015mJmol-'^K-^T-^. By taking ^ = 2ll, we obtain 
the following values: a — 2AA9mJmol^^K~^ and 62 — 
2.239mJmor'^K-^. It is known that a = 2.52jn/Tc, we 
further obtain 7„ — 26.63mJmol~^K~^ . In addition, it 
is predicted that bi ~ k'jn/VHc2 which yields also 7„ = 
23.3m Jmol~^K~^ if taking k = 1 and Hc2 = 4'o/'^£.'^- 
Two different approaches lead to rather close values of 
7„. This value is also close to that found in y — 1230. 

The nice fitting in Fig. 3 with only the first term of 
eq.5 suggests that the core region has very small contri- 
bution to the DOS since otherwise the last term jcoreH 
should be sizeable. This implies that the low energy DOS 
inside the vortex core is very small, ruling out the pres- 
ence of ZBCP as expected by the BdG theory for a d- 
wave superconductor 0, 0|. This conclusion is further 
strengthened in the following by the fact that the core 
size effect has a sizeable influence on the DOS due to 
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FIG. 3: DOS induced by Doppler shift at zero K. The solid 
line is a theoretical curve A7 = 1.9VH{1 - 0.0067//). The 
remarkable consistency between data and theory based on 
d^2_y2 symmetry indicates that we don't have any second 
component such as id^y or is at fields ranging from to 12 T. 
The inset shows the slope d'y/dT of the straight lines shown 
in Fig. 2 in low temperature region. The dashed line is a linear 
fit to the data at fields above one tesla. The intercept and 
the slope of the dashed line give rise to the pre-factors of the 
second and third terms in eq.5 



Doppler shift, but the contribution from inner part of 
the core is small. In order to further test our idea ( eq.5 
), we present the scaling of A7r— (62 — a)?"^ vs. T\/H in 
Fig. 4 with 62 — a derived above. The quality of our scal- 
ing is remarkable. This can be easily understood from 
eq.5 when the correction 7rf^i//0o(= 0.08 << 1 at 12 T 
) to the first term is very small. The slight scattering in 
the high temperature and field region is due to the noise 
of the data. Worthy of noting is that to have this nice 
scaling we need to take ^core ~ 0; again showing a small 
contribution from the inner vortex core. The solid line 
is a theoretical curve from eq.5 by using the parameters 
derived above. Both the nice scaling and the consistency 
between the experimental data and cq.5 suggest that the 
Simon-Lee scaling law can be reconciled by considering 
the vortex core size effect. Meanwhile the scaling accord- 
ing to the original Simon-Lee theory is also presented in 
the inset of Fig. 4. showing a rather poor scaling quality. 

In low temperature region, our analysis indicates that 
the field induced DOS can be well described by Volovik's 
theory or Simon-Lee scaling law although a correction 
due to the core size effect is needed. This means that the 
prerequisite for the theory, i.e., the dj.2_y2 pairing sym- 
metry is well satisfied. Therefore it naturally rules out 
the possible presence of a second order parameter like 
idxy or is either due to overdopingj^ or due to the field 
effect jl^ in the present overdoped sample. Meanwhile 
another seemingly contradicted phenomenon is that the 
vortex core region contributes very little ( at least much 
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FIG. 4: Scaling of the data Cscai = C{H) - C{H = 0) + (a - 
62) vs. Ty/H. The solid line is a theoretical fit according 
to eq.5 with fitting parameters derived in the text. A perfect 
scaling and a nice consistency between data and the scaling 
law with core-size correction is clear. The inset shows the 
scaling according to Simon and Lee ( eq.l ) where Cscai = 
C{H)-C{H = 0)+aT^. 



smaller than that induced by the Doppler shift if the 
super-current would flow in the same area ) to the to- 
tal DOS. This suggests that there is no ZBCP within 
the normal core since otherwise the normal core region 
should give rise to a sizable signal. Our conclusion is con- 
sistent with the tunnelling results[ll[ll[ll|23,|23 and 
certainly clears up the concerns about the surface con- 
ditions in the STM measurement since our specific heat 
data reflect the bulk property. Recent results from NMR 
also show the absence of ZBCP inside the vortex core|2^. 
In this sense our data together with the earlier NMR data 
present a bulk evidence for an anomalous vortex core. 
Interestingly it is widely perceived that the normal state 
in overdoped region shows a Fermi liquid behavior even 
when the superconductivity is completely suppressed jl^. 
Clearly the mean-field frame of BdG theory based on the 
conventional d-wave superconductivity is not enough to 
interpret the anomalous vortex core state. 

In conclusion, the quasiparticle density of states ( DOS 
) due to Doppler shift in the mixed state of an overdoped 
La2^xSrj.Cu04 ( x = 0.22 ) single crystal is found to 
deviate from the proposed scaling law of Simon and Lee. 
However, this law is reconciled if one considers the core 
size effect . The contribution from the inner vortex core 
is small comparing to that due to the Doppler shift in the 
same area. Our results suggest the absence of the ZBCP 
in the vortex core although it is expected by the Bogoli- 
ubov de-Gennes theory for a d-wave superconductor. 
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